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Using 180 pb _1 of data collected with the CDF II detector at the Tevatron, we measure the first 
two moments of the hadronic invariant mass-squared distribution in charmed semileptonic B decays. 
From these we determine the non-perturbative Heavy Quark Effective Theory parameters A and 
Ai used to relate the B meson semileptonic branching ratio to the CKM matrix element \V c b\- 
For a minimum lepton momentum of 0.7 GeV/c in the B rest frame we measure the first two 
moments of the D** — > D^n component to be (m 2 D **) = (5.83 ± 0.16 sta t ± 0.08 syst ) GeV 2 /c 4 and 
((m 2 D „ - (m 2 D *,)) 2 ) = (1.30 ± 0.69 sta t ± 0.22 syst ) GeV 4 /c 8 . Combining these with the discrete 
mass terms from the D and D* mesons, we find the total moments to be (M| c ) — fri 2 D = (0.467 ± 
0.038 stat ±0.068 syst ) GeV 2 /c 4 and {{M 2 Xc - (M 2 Xc )) 2 ) = (1.05 ± 0.26 sta t ± 0.13 syst ) GeV 4 /c 8 , where 
TUd is the spin-averaged D mass. The systematic error is dominated by the uncertainties in the 
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world-average branching ratios used to combine the D, D* , and D** contributions. The analysis 
makes no assumptions about the shape or resonant structure of the D** D^n invariant mass 
distribution. 

PACS numbers: 13.20.He, 12.15.Hh, 12.39.Hg 



I. INTRODUCTION 



In order to constrain the length of the side opposite 
the angle (3 in the Cabibbo-Kobayashi-Maskawa (CKM) 
unitarity triangle, a precise measurement of the ratio 
I | / 1 V^b | is needed. The matrix element \V c b\ is gen- 
erally extracted from semileptonic B decays. Currently, 
the most precise method is based on the measurement 
of the inclusive semileptonic partial width into charm, 
r s ; = T(B — ► X c lvi). The Operator Product Expan- 
sion (OPE) applied to Heavy Quark Effective Theory 
(HQET) relates the experimental determination of T s ; to 
I Vcb CH3 The relationship takes the form of an expan- 
sion in inverse powers of the B mass, tub- At each order 
in the expansion, new free non-perturbative parameters 
enter: one (A) at order 1/ms, two (Ai and A2) at order 
l/m B , six at order l/m B , etc. In order to extract \V c b\ 
from T s i some external information on these parameters 
is needed. 

The same theoretical framework that predicts the value 
of r s ; predicts the value of any weighted integral of 
the differential rate dT s i/dsH, provided the weight is a 
smooth function of sh = M| . Using weight functions 
(sh — m 2 D ) and (s# — (sh)) 2 , with Sj = 0.25m/) + 
0.75m/)*, one can define the first two moments of the 
hadronic mass distribution: 



cntial mass-squared spectrum can be written as: 
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where T s i is now the inclusive B~ semileptonic width, 
To and T* are the exclusive B~ partial widths to D°l~Vi 
and D*°l~vi respectively, and f**(sn) is the normal- 
ized hadronic invariant mass-squared distribution in the 
D**° channel. We use world-average values of Tq/T^j. 
T*/T s i, to £)0 and 777£>*o from the Particle Data Group [3j 
and concentrate on measuring f**(su)- In this way, we 
have only to measure the invariant mass distribution for 
the D**° component without having to determine the 
D , D*° components or the relative normalizations be- 
tween those and the D**° channel. The D**° spectrum 
is not well known, and includes, at least, two narrow and 
two wide states, together with a possible non-resonant 
D^*'mr contribution. The measurement of the D**° spec- 
trum is the main task of this analysis. We assume that 
the D^ttI~vi decays of B~ saturate the difference be- 
tween its inclusive semileptonic decay rate and the sum 
of its exclusive decay rates to D°l~Vi and D*°l~Vi. We 
neglect all modes with additional pions in the final state, 

as well as D**° -> D { & * )+ K~ . 

Only Z?(*) + 7T~ decays (charge conjugated channels are 
implicitly included throughout the paper) are recon- 
structed. Contributions to the sh distribution from 
decays with neutral particles are included by applying 
isospin factors to the charged modes. Feed-down from 
one channel to another due to unmeasured neutral parti- 
cles is subtracted statistically using the data themselves 
and isospin relations, as explained in Section ITTT1 



which are simply the shifted mean and variance of the 
M\ distribution in semileptonic charmed decays of B 
mesons. The moments are not sensitive to |Vc&|, but 
they are more sensitive to the non-perturbative parame- 
ters of HQET than r s ; itself is. Therefore, measuring the 
moments provides a useful constraint on the HQET pa- 
rameters which improves the overall precision on | V c b \ as 
determined from r s ;. This is the purpose of this analysis. 
Since A2 is well determined from the values of the hyper- 
fine mass splittings in the B and D meson systems 0, 
only A and Ai are studied here. 

The sh distribution in B~ — > X®l~Vi decays can be 
split into three contributions corresponding to X® = 
D°, D*°, D**°. Here D**° stands for any neutral charmed 
state, resonant or not, other than D°, D*°. The diffcr- 



II. DATA ANALYSIS 

The analysis uses a data sample of pp collisions at 
= 1.96 TeV with an integrated luminosity of about 
180 pb" 1 , collected between February 2002 and August 
2003 with the upgraded Collider Detector at the Fer- 
milab Tevatron (CDF II). A description of the detector 
can be found in [4j. The relevant components for this 
analysis include a tracking system composed of a silicon 
strip vertex detector (SVX II) surrounded by an open 
cell drift chamber system (COT). The SVX II detector 
comprises five concentric layers of double-sided sensors 
located at radii between 2.5 and 10.6 cm, while the COT 
provides 96 measurements (including axial and stereo) 
out to a radius of 132 cm. The central tracking sys- 
tem is immersed in a 1.4 T solenoidal magnetic field. 
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Two sampling calorimeters surround the magnetic coil. 
A set of proportional chambers inside the electromagnetic 
calorimeter provides information on the shower profile for 
use in electron identification. Muon candidates are iden- 
tified with two sets of multi-layer drift chambers, one 
located outside the calorimeters and the other behind an 
additional 60 cm-thick iron shield. 

Decays B — > D^ + n^l^X : where I stands for elec- 
tron or muon, were recorded using a trigger that re- 
quires a lepton I and a track displaced from the inter- 
action point 0- The lepton and the displaced track 
must have transverse momentum pt in excess of 4 GeV/c 
and 2 GeV/c respectively. The displaced track's im- 
pact parameter with respect to the beamline has to ex- 
ceed 120 /xm and be below 1 mm. Events which pass 
the trigger are recorded for further analysis. Only well- 
reconstructed tracks with px > 0.4 GeV/c are retained. 
Track parameters are corrected for the ionization energy 
loss appropriate to the mass hypothesis under considera- 
tion. A Monte Carlo sample of B — > D**lui events based 
on the ISGW2 [|J and Goity-Roberts matrix elements 
and including a detailed simulation of the CDF II de- 
tector based on the GEANT 8] package has been used 
throughout the analysis. In accordance with our assump- 
tion, only D** — > DWtt decays are generated. 

Events with D* + (^ D a n+)l~ and D + l~ combina- 
tions are reconstructed in the decay channels D° — > 
K~ir + , K~tt + it~it + , K~ir + ir° and D + — > K~tt + it + . 
Tracks with the appropriate charge combination are re- 
quired to be consistent with a common vertex in three 
dimensions. One of the tracks in the vertex must fulfill 
the displaced trigger requirements. Suitable ranges are 
selected in the D° (1.84-1.89 GeV/c 2 ) and D+ (1.84- 
1.89 GeV/c 2 ) mass distributions. For the D* + chan- 
nel, an additional charged track (7r+) is required, such 
that the M(D°ir+) — M(D°) mass difference lies between 
0.142 and 0.147 GeV/c 2 . The D° -> R-tt+tt chan- 
nel is reconstructed from the satellite peak in the K~w + 
mass distribution (1.50-1.70 GeV/c 2 ). In this case the 
M(K~ir + ir+) — M(K~ir + ) mass difference is required 
to be between 0.142 and 0.155 GeV/c 2 . Duplicate re- 
moval is performed in the D° — > K~ir + ir~ir + channel: 
when two D* + candidates share all five tracks, differing 
only in the kaon mass assignment, the candidate with 
the K~7t + it~tt + mass closer to the nominal D° mass is 
retained. No attempt has been made to further identify 
kaons and pions. 

After the selection of D* + l~ and D + l~ combina- 
tions, we obtain 3890 ± 63, 2994 ± 57, 6638 ± 98 and 
14416 ± 202 signal events in the K~ir + , K~ir + ir~ir + , 
K~ir + ir° and K~ir + ir + channels, respectively. Combi- 
natorial backgrounds, estimated from sidebands of the 
M (D+) and M (D°n+)- M(D°) distributions, have been 
subtracted. The quoted yield in the D + channel has 
been rescaled by a factor 0.96 to account for the back- 
ground from Df — > K + K~tt + decays where the K + is 
assigned the pion mass. Figure ^ shows on the left the 
M(D°tt+) - M(D°) distributions for D° decaying into 



either K~?r + or i-C ~7r + 7r~7r + , and for D a — > _ftT~7r + 7r , 
while the D + — > A"~7r + 7r + mass distribution is plotted 
on the right. 

The D^ + l~ vertex (the B vertex) is reconstructed 
in three dimensions and required to be at least 500 /im 
away from the beam line. An additional pion (tt^) is 
then added to create full D + ^*\~l~ candidates. The 
Tii's trajectory is required to be at most 2.5 standard 
deviations away from the B vertex, and at least three 
standard deviations away from the beam line. These cuts 
were optimized using the B — > D**lvi Monte Carlo for 
the signal, and wrong-sign 7rJ<J _ combinations in data 
for the background. The measured mass distributions in 
the Z?* + 7r»» and D +r K stst channels are shown in Fig. [21 



III. BACKGROUND AND EFFICIENCY 
CORRECTIONS 

The most important background sources are combi- 
natorial background under the D + and D* mass peaks, 
prompt tracks (from fragmentation or the underlying 
event) that fake 7r„» candidates, and feed-down from 
D* + — > D + tt° decays into the D + channel. Data side- 
bands are used to assess combinatorial background un- 
der the D + mass and D* + -D° mass-difference peaks, 
and wrong-sign tt^I" combinations in data to character- 
ize the prompt background to the 7r»» candidates. The 
wrong-sign pion-lepton sample is subtracted from the 
right-sign sample, after performing side-band subtraction 
in both. A possible difference between the rate of prompt 
background with the right and the wrong 7r„* charge has 
been studied with a sample of fully reconstructed B de- 
cays (B- -> J/VK-, B° -> L>M-7T + , B~ -> D°tt~) 
and found to be at most 4%. This has been included in 
the systematic error. The bias in the background sub- 
traction introduced by using the same wrong-sign sam- 
ple for both the optimization of the selection and the fi- 
nal background subtraction has been studied using boot- 
strap 9j copies of the data and found to be smaller than 
15% of the statistical error on the moments. This up- 
per bound on the bias has been added as an additional 
systematical error. 

Since we do not reconstruct neutral particles, events 
with a decay B~ — > D* + tt^I^Vi with the D* + decay- 
ing into D + tt° constitute an irreducible background to 
the signal channel B~ — > D + tt^I^Vi. Using isospin 
symmetry, the background rate and D + n^ lf invariant 
mass can be obtained from the D ^^ invariant mass in 
B~ -> D*+n-J-Vi decays with D*+ -> D°tt+, which 
we fully reconstruct, after correcting for the relative ef- 
ficiency using Monte Carlo. A small physics background 
(around 1% in rate), coming mostly from B — > DgD^ 

decays with the fi* decaying semileptonically, is sub- 
tracted using Monte Carlo predictions. Background from 
tracks faking a lepton has been studied by looking at 
wrong-sign D^ + l + combinations. It has been found that 
the background subtraction procedure outlined above 
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FIG. 1: Left: Mass difference Am* = M{D°ir+) - Af(D°) in the D° -> K~ty + and D° -* R-ty+ty+tt- channels (narrow 
peak) and Am* = M(K~ Tv + irf) — M(K~TT ) in the D° -> 7r + 7r channel (broad peak). Right: Mass distribution for the 
D + — » A' _ 7r + 7r + channel. The signal areas are shaded; sideband regions are also indicated. 
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FIG. 2: Sideband-subtracted invariant mass distribution for the D* + ir,, channels (left) and for the D + ir,, channel (right). 
The mass regions shown are limited at 3.5 GeV/c 2 for illustration only. No explicit mass cut is applied in the analysis. 



effectively removes any such background. Finally, the 
background from B — > D^-ktv t decays has been stud- 
ied with Monte Carlo and found to have a negligible effect 
on the determination of the moments. 

Since only the shape of the mass-squared distribution 
for the D** component (/*'(s H ) in Eq. (0)) is being mea- 
sured, the relevant efficiency corrections are those that 
can bias the mass-squared distribution, along with the 
relative efficiency for the D* + and D + components of 
the D** piece. Both efficiencies are obtained, as a func- 
tion of the mass M(D**), from the Monte Carlo simula- 
tion. We have checked the Monte Carlo relative efficiency 



predictions as a function of the transverse momen- 
tum by applying the selection cuts to decay tracks 
from D and D* mesons. The efficiency variation in the 
data agrees well with the simulation. Where small dif- 
ferences have been found, corrections have been derived 
from these Monte Carlo-data comparisons. 

To compare the moments with theoretical predictions, 
they must be measured with a well defined cut on p*, the 
lepton momentum in the B rest frame. Since we do not 
attempt to measure the boost of the B, we cannot access 
p* directly in data. Instead, acceptance corrections are 
derived from Monte Carlo that turn our gradual trigger 
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turn-on as a function of p\ into a sharp threshold at p\ = 
0.7 GeV/c, thereby correcting our measurement of the 
moments to a cut p* > 0.7 GeV/c. The value 0.7 GeV/c 
was chosen in order to minimize the correction. Because 
of the negative correlation between lepton momentum 
in the B rest frame and D** mass, the correction itself 
can depend on the detailed D** mass spectrum in Monte 
Carlo. In order to assess the possible systematic error, 
the default B decay model has been compared to a naive 
phase-space B scmilcptonic decay model. The differences 
in the ensuing correction factors as a function of M(D**) 
are considered as systematic errors. 



IV. RESULTS 

The Z)(*) + 7r _ mass distribution is shown in Fig. [3] af- 
ter background subtraction and efficiency and acceptance 
corrections. The first and second moments of the D** 
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FIG. 3: Fully corrected invariant mass distribution 
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background subtracted and corrected for mass-dependent and 
D* J D + relative efficiency corrections. The plotted errors take 
into account all corrections and subtractions. 



component of the mass-squared distribution, m\ and mi , 
are determined by simply computing the mean and vari- 
ance of the distribution shown in Fig. |3 without any 
assumption about the shape or rate of its several compo- 
nents: 



mi 
m 2 



(5.83 ± 0.16 sta t ± 0.08 syst ) GeV 2 /c 4 

(ml, t )f) = 

(1.30 ± 0.69 sta t ± 0.22 syst ) GeV 4 /c 8 , 



with a 61% positive correlation. The full moments of 
the hadronic mass-squared distribution, Mi and Mi, are 
determined by combining mi and mi with the D and D* 



T 2 \ — 2 
m n 



2 1 A 



pieces, obtained from world-average values 0: 
Mi = (M|J 

(0.467 ± 0.038 stat ± 0.019 oxp ± 0.065 BR ) GeV 2 / 
Mi ee ((M 2 Xr - (MiJ) 2 ) = 

(1.05 ± 0.26 stat ± 0.08 cxp ± O.IObr) GeV 4 /c 8 , 

with a 69% positive correlation between Mi and Mi. 
Here "BR" refers to the uncertainty coming from the 
branching ratios needed for the combination of the D, 
D* and D** pieces. For the exclusive branching ratios 
to D and D* , all available information |3( coming from 
charged and neutral B decays has been combined using 
isospin invariance, leading to Fo/r s ; = 0.203 ±0.015 and 
r»/r s; = 0.550 ± 0.026 with about 30% positive correla- 
tion. The isospin-related partial widths (not the branch- 
ing ratios) of B + and B° are assumed to be identical. 

Finally, using the predictions in 2] , the HQET parame- 
ters A and Ai are determined. After applying constraints 
on the other HQET parameters coming from the known 
B and D hyperfine mass splittings we find, in the pole 
scheme: A = A polc = (0.397 ± 0.078 stat ± 0.027 cxp ± 



0.064 BR ± 0.058 th co) GeV and \\ 



olc 



(-0.184 ± 



0.057 sta t ± 0.017 oxp ± 0.022 BR ± 0.077 thc o) GeV 2 , with 
a 79% negative correlation. Similarly we extract the 
equivalent HQET parameters in the IS scheme: 



M T /2 - A 16 = (4.654 ± 0.078 sta t ± 0.027 cxp ± 0.064 BR ± 
0.089 thco ) GeV and A} s = (-0.277±0.049 stat ±0.017 cxp ± 
0.022 BR ± 0.094 thc o) GeV 2 , with a 77% positive correla- 
tion. 

The statistical and systematic errors in the extraction 
of the D** moments, the full moments, and the HQET 
parameters in the pole-mass scheme are presented in Ta- 
ble n Statistical errors dominate the measurements of 
mi and mi while experimental systematic errors are all 
smaller. The main experimental systematics are com- 
puted from the differences in the results when we apply 
or omit a correction for the ~ 60 MeV/c 2 mass resolution 
in the K~7r + n° channel or the M(Z?**)-dependent effi- 
ciency correction from data. Similarly, we considered the 
difference in results obtained using Monte Carlo efficiency 
and acceptance corrections calculated with ISGW2 and 
Goity-Roberts matrix elements or phase space as a sys- 
tematic error. Other experimental systematics include 
uncertainties on the level of the prompt background 
(studied with a fully reconstructed B sample), a possible 
bias due to having used the same data sample to model 
the background in the optimization process and to sub- 
tract the background from the data (studied by repeat- 
ing the optimization on bootstrap copies of the data), 
and uncertainties on D branching fractions used in the 
analysis. 

Uncertainties in the inclusive and exclusive semilcp- 
tonic B branching ratios become important when com- 
bining mi and mi with the D a and D*° pieces to obtain 
Mi and Mi, the moments of the entire charm mass dis- 
tribution. Theoretical uncertainties become dominant in 
the extraction of the HQET parameters A and Aj.. The 
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Error 


Ami 

(GeV 2 /c 4 ) 


Am2 
(GeV 4 /c 8 ) 


AMi 

(GeV 2 /c 4 ) 


AM 2 
(GeV 4 /c 8 ) 


AA 

(GeV) 


AAi 
(GeV 2 ) 


Statistical 
Total systematic 


0.16 
0.08 


0.69 
0.22 


0.038 
0.068 


0.26 
0.13 


0.078 
0.091 


0.057 
0.082 


Mass resolution 


0.02 


0.13 


0.005 


0.04 


0.012 


0.009 


Efficiency (data) 


0.03 


0.13 


0.006 


0.05 


0.014 


0.011 


Efficiency and acceptance (MCI 


0.06 


0.05 


0.016 


0.03 


0.017 


0.006 


Background scale 


0.01 


0.03 


0.002 


0.01 


0.003 


0.002 


Background bias 


0.02 


0.10 


0.004 


0.03 


0.006 


0.006 


Physics background 


0.01 


0.02 


0.002 


0.01 


0.004 


0.002 


D + 1 D* + branching ratios 


0.01 


0.02 


0.002 


0.01 


0.004 


0.002 


D + /D* + efficiency 


0.02 


0.03 


0.004 


0.01 


0.005 


0.002 


B semileptonic branching ratios 






0.065 


0.10 


0.064 


0.022 


Pi 










0.041 


0.069 


Ti 










0.032 


0.031 


a a 










0.018 


0.007 


m b ,m c 










0.001 


0.008 


Choice of p ; * cut 










0.019 


0.009 



TABLE I: Statistical and systematic uncertainties in the measurements of the D** and full moments and in the extraction of 
A and Ai in the pole scheme. 



largest contribution to the theoretical systematic error 
is that estimated by varying the unknown third order 
HQET parameters in the ranges p\ = ^(0.5 GeV) 3 ± 
|(0.5 GeV) 3 , Ti = (0.0 GeV) 3 ± (0.5 GeV) 3 . Finally, 
acceptance corrections have been computed for two al- 
ternative Pi cuts, 0.5 and 0.9 GeV/c. The moments ob- 
tained this way are different physical quantities and nu- 
merically different from those obtained for the default p*{ 
cut at 0.7 GeV/c. However, if HQET describes the data, 
they should all lead to compatible values of the HQET 
parameters. The three sets of parameters are found to 
be equal within errors. Differences between them have 
been considered as additional systematic uncertainties. 

In summary, we have presented a measurement of the 
first two moments of the hadronic mass-squared distri- 
bution in semileptonic B decays to charm by combin- 
ing our measurement of the D^ + tt^ mass spectrum 
with the known masses and branching ratios to Dlv and 
D*lv taken from the Particle Data Group compilation Q. 
These channels together are assumed to fully account 
for the inclusive semileptonic decay width of B mesons 
to charm. The moments are then used to extract the 
two leading Heavy Quark Effective Theory parameters, 
A and Ai , in both the pole mass and the IS mass schemes. 
Within HQET the results are in agreement with previous 
determinations at e + e~ machines [Toj| . mostly at higher 
values of and our precision in the moments is com- 
parable or slightly better. The experimental techniques 
and systematic uncertainties are very different from those 
at electron-positron colliders. Our measurements, when 
combined with moments of the lepton energy in semilep- 



tonic B decays [Tl| and photon energy in b — > 57 transi- 
tions |l2j measured elsewhere, will help pinpoint the val- 
ues of the non-perturbative HQET parameters • This 
will sharpen the experimental determination of \V c b\ and 
allow tests of the underlying assumptions of the theoret- 
ical framework, such as quark-hadron duality. 



ACKNOWLEDGMENTS 

We thank the Fermilab staff and the technical staffs of 
the participating institutions for their vital contributions. 
It is a pleasure to thank Zoltan Ligeti for several en- 
lightening discussions regarding the HQET predictions. 
This work was supported by the U.S. Department of En- 
ergy and National Science Foundation; the Italian Isti- 
tuto Nazionale di Fisica Nucleare; the Ministry of Educa- 
tion, Culture, Sports, Science and Technology of Japan; 
the Natural Sciences and Engineering Research Council 
of Canada; the National Science Council of the Repub- 
lic of China; the Swiss National Science Foundation; the 
A. P. Sloan Foundation; the Bundesministerium fur Bil- 
dung und Forschung, Germany; the Korean Science and 
Engineering Foundation and the Korean Research Foun- 
dation; the Particle Physics and Astronomy Research 
Council and the Royal Society, UK; the Russian Foun- 
dation for Basic Research; the Comision Interministerial 
de Ciencia y Tecnologi'a, Spain; and in part by the Euro- 
pean Community's Human Potential Programme under 
contract HPRN-CT-2002-00292, Probe for New Physics. 



[1] A. Falk, M. Luke, M. Savage, Phys. Rev. D 53, 2491 (1997); A. Falk, M. Luke, Phys. Rev. D 57, 424 (1998). 

(1996); M. Gremm, A. Kaputsin, Phys. Rev. D55, 6924 



9 



[2] C. W. Bauer, Z. Ligeti, M. Luke, A. V. Manohar, Phys. 

Rev. D 67, 054012 (2003). 
[3] S. Eidelman et al, Phys. Lett. B592, 1 (2004). 
[4] D. Acosta et al. (CDF Collab.), Phys. Rev. D68, 072004 

(2003) . 

[5] B. Ashmanskas et al, Nucl. Instrum. Moth. A518, 532 

(2004) . 

[6] D. Scora, N. Isgur, Phys. Rev. D 52, 2783 (1995). 

[7] J. L. Goity, W. Roberts, Phys. Rev. D 51, 3459 (1995). 

[8] R. Brun, F. Carminati, CERN Programming Library 
Long Write-up W5013 (1993). 

[9] B. Efron, Ann. Stat. 7, 1 (1979); A. Stuart, J. K. Ord, 
Kendall's Advanced Theory of Statistics, Vol. 1, 5th Ed., 
(Oxford University Press, New York, 1987). 



[10] D. Cronin-Hennessy et al. (CLEO Collab.), Phys. Rev. 

Lett. 87 251808 (2001); M. Battaglia et al, Phys. Lett. 

B556, 41 (2003); B. Aubert et al. (BABAR Collab.), 

Phys. Rev. D 69, 111103 (2004). 
[11] H. Mahmood et al, (CLEO Collab.), Phys. Rev. D 67, 

072001 (2003); M. Battaglia et al, Phys. Lett. B556, 41 

(2003); B. Aubert et al. (BABAR Collab.), Phys. Rev. D 

69, 111104 (2004). 
[12] S. Chen et al (CLEO Collab.), Phys. Rev. Lett. 87, 

251807 (2001). 

[13] C. W. Bauer, Z. Ligeti, M. Luke, A. V. Manohar, 
M. Trott, Phys. Rev. D 70, 094017 (2004). 



